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Abstract 

The Method of Invariant Grid (MIG) is a model reduction technique based on the concept of slow invariant manifold 
(SIM), which approximates the SIM by a set of nodes in the concentration space (invariant grid). In the present 
work, the MIG is applied to a realistic combustion system: An adiabatic constant volume reactor with i?2-air at 
stoichiometric proportions. By considering the thermodynamic Lyapunov function of the detailed kinetic system, 
the notion of the quasi-equilibrium manifold (QEM) is adopted as an initial approximation to the SIM. One- and 
two-dimensional discrete approximations of the QEM (quasi-equilibrium grids) are constructed and refined via the 
MIG to obtain the corresponding invariant grids. The invariant grids are tabulated and used to integrate the reduced 
system. Excellent agreements between the reduced and detailed kinetics is demonstrated. 
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1. Introduction 

Accurate modeling of reactive flows of hydrocar- 
bon fuels requires the solution of a large number of 
conservation equations as dictated by detailed reac- 
tion mechanism. In addition to the sometimes pro- 
hibitively large number of variables introduced, the 
numerical solution of the governing equations has to 
face the stiffness introduced by the fast time scales 
of the kinetic terms. These issues make computa- 
tions of even simple flames time consuming. On the 
other hand, the dynamics of complex reactive sys- 
tems is often characterized by short initial transients 
during which the solution trajectories approach low- 
dimensional manifolds in the concentration space, 
known as the slow invariant manifolds (SIM). The 
remaining dynamics lasts much longer and evolves 
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along the SIM towards the equilibrium state. The 
construction of the SIM allows to establish a simpli- 
fied description of a complex system by extracting 
only the slow dynamics and neglecting the fast. As 
a result, the detailed large set of equations can be 
reduced to a much smaller system without a signifi- 
cant loss of accuracy. For this reason, much effort is 
devoted to develop automated model reduction pro- 
cedures based on the concept of SIM. The methods 
of Intrinsic Low Dimensional Manifold (ILDM) pQ 
and Computational Singular Perturbation (CSP) [5] 
represent the two examples of this family of meth- 
ods. 

In the sequel, the Method of Invariant Grid (MIG) 
introduced in \3\ as a computational realization of 
the Method of Invariant Manifold (MIM) g] is used 
for the first time to study a combustion problem. 
The Quasi-Equilibrium Grid algorithm introduced 
in [5] is adopted to obtain a first approximation of 
the SIM, which is afterwards refined via MIG it- 
erations. The data delivered by that procedure are 
stored in tables and used to integrate a smaller and a 
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less stiff reduced system. The results obtained with 
the reduced system are compared with the detailed 
one, and excellent agreement is found for the dy- 
namics of all species, including the radicals, and of 
the temperature. 



the vector state ip is not explicitly known and, in 
general, the evaluation of the right-hand side of ([3]) 
requires implicit solution of the second equation in 
©■ 



2. Problem setup 

Consider a reactive system where n chemical 
species x\,...,x n participate in a complex mech- 
anism with r reactions. Let a generic reversible 
reaction step be written as: 

a„lX\ + . . . + a sn X n ^ f3 B lXl + . . . + PsnXn, (1) 

where s = 1, . . . , r is the reaction index and the inte- 
gers a S i and j3 S i are the stoichiometric coefficients of 
the reactants and products in the reaction s, respec- 
tively. Let the corresponding stoichiometric vectors 
be ac s = (a sl , a sn ), (3 S = {(3 sl , (3 sn ) and 7 S = 
(3 S — a. s . The temporal evolution of a homogeneous 
reactive system is described by a system of ordinary 
differential equations (ODEs) 

dtp 
~dt 

where tp is a set of variables that characterizes the 
thermo-chemical state of the system. In the sequel, 
we consider an adiabatic constant volume reactor 
where the density and the mixture-averaged spe- 
cific energy are fixed. Let [X,] denote the molar 
concentration of species i. Once the vector xp = 
(p, e, [Xi], [X„]) T is taken to describe any state 
of the system, then @ reads: 



(2) 



g (i/0 = 0, 0, ]T 7.(1)0-, ]T ls (n)Q s 



dip 
~dt 

\ s=l s=l / 

(3) 

where the superscript T denotes the transposition. 
The rate of reaction s, and the mean internal 
energy can be expressed as 

n s = Qf - 0" 

n n 

^ = k (t) n , nj = k (t) n [x t f 

i=l i=l 

n 
i=l 

(4) 

where fc+ ,k s ,ei,Yi are the forward and inverse re- 
action rate constants of the reaction s, and the spe- 
cific internal energy and the mass fraction of species 
i, respectively. The dependence of temperature T on 



2.1. Thermodynamic Lyapunov function 

Here we assume that the kinetic system ([2]) de- 
scribes the evolution of a chemical system towards a 
unique equilibrium state. Moreover, we assume that 
there exists a strictly convex function dependent on 
the state vector tp that decreases monotonically in 
time under the dynamics of the system ([2]). Such a 
function G is a global Lyapunov function of the sys- 
tem, and it reaches the global minimum at the equi- 
librium point. 

In a closed system under constant energy e and 
density p, the value of specific mixture-averaged en- 
tropy s must increase monotonically starting from 
any initial condition. Therefore, the function G — 
—s 



G = 



E 



Sl {T)-R\n{Xi)-R\n(j^) 



Xi 



w 



_(5) 

is a Lyapunov function of system ([3]), where W is 
the mean molecular weight, Sj and Xi are the en- 
tropy and the mole fraction of species i, respectively, 
R is the universal gas constant while p and p re f 
are the mixture total pressure and a given reference 
pressure, respectively. Assuming that d elements are 
involved in the reaction, we can construct another 
Lyapunov function 

d / n \ n 

G = G+J2 A^ftW ]+xY,W i [X i ], (6) 



k=l 



i=l 



where pki represent the number of atoms of the k-th 
chemical element in species i and Wi is the molecu- 
lar weight of species i. The fixed parameters Xk and 
A are chosen such that VG| g = at the equilib- 
rium, with the gradient of G computed under fixed 
e. Because of the conservation of atoms and density, 
the time derivative of © is non-positive: 



dG 
~dt 



dG 
~dt 



<0, 



dNk 
dt 



0. 



dp 
~dt 



0. 



Here, Nk — E"=i Mfci [Xi] denotes the total number 
of atoms k in the system and p — E™=i ^ ■ 
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3. Theoretical background 

In the present work, we consider the method of 
invariant grid (MIG) for combustion applications. 
In the sequel, we refer to a discrete collection of 
points in the concentration space as a grid. Start- 
ing from an initial approximation, the MIG approx- 
imates the slow invariant manifold (SIM) of system 
© by means of the invariant grid, and refines the 
initial grid iteratively. At each iteration, the correc- 
tions are obtained by solving a system of linear equa- 
tions at every node. The set of refined nodes forms 
a new grid, which approximates the SIM more accu- 
rately than the initial grid. Next iterations are then 
carried out until a termination criterion is satisfied. 
Finally, the reduced dynamics on the invariant grid 
is obtained with the help of a proper parameteriza- 
tion. 

More details about the theoretical background of 
the MIG and its implementation can be found in 
Refs. [3I4I5I6I7] . In the sequel, we focus on the as- 
pects related to MIG application to non-isothermal 
cases. 

3.1. MIG algorithm 

Let us consider a g-dimensional manifold f2 in 
the n-dimensional concentration space. We assume 
that q variables £ l are associated with each point 
of f2 by a smooth function F{£} , ■ ■■,t; q ) which maps 
the variables £ l into the concentration space. The 
tangent plane r to the manifold at a given point of 
fi is spanned by q vectors: 



OF/dC 



1, 



,q. 



(7) 



Any n-dimensional vector x can be projected onto 
the tangent plane r by introducing a projector P, 
such that the projected n-dimensional vector P(x) 
belongs to r and P (P (x)) = P (x). The manifold 
fl is invariant with respect to ([3]) if any solution 
trajectory starting on fl proceeds towards the equi- 
librium state along this manifold. In other words, 
the difference between the vector field g and its pro- 
jection on the tangent space defines the invariance 
condition 

(I P)g = (8) 

which must be satisfied at every point of f2, with I 
denoting the identity matrix. 

MIG is an iterative procedure which aims at re- 
fining an initial set of points in the concentration 
space (initial grid) and delivering the invariant grid 



that describes the slow invariant manifold. More 
generally, MIG regards the invariance condition f8} 
as an equation and solves it by Newton-like itera- 
tions. Consider a set of nodes Q in the concentration 
space (grid) approximating the g-dimensional SIM. 
Assume that at any node of Q the q tangent vec- 
tors Ui can be approximated by finite differences, 
and the construction of the projector P is defined 
(a specification of P will be given below). The grid 
Q is considered as invariant when the left-hand side 
of ([8|) is sufficiently small with respect to the vector 
field g at each node. 

In order to take into account the conservation of 
elements and density, we consider the (d + 1) x n 
matrix: 



M 



Mil ■ ■ • Mir 



Mrfi 



(9) 



Wi ... W n 

Let the vector set {bi , . . . , bh} be selected as a basis 
of the /i-dimensional intersection of the null space of 
the projector P and the null space of M. Where h 
is the dimension of the latter intersection subspace. 
Following the MIG approach, once an initial grid Go 
(in general, non-invariant) is given, it can be refined 
by solving at each node the following system of h 
equations |6j 



h 

£« 

i=l 



bjJbt 



(jbf)] = [P(f)bj-fb 



(10) 

for the unknowns 8i. The vector field / represents 
the time derivatives of species concentrations, / = 
(d[Xi]/dt), and J = [df/d[Xi]} denotes the first 
derivative matrix (Jacobian) of /. Once the alge- 
braic system (fT0|) is solved at each point X° = 
[X®]) of Go, a new point X 1 is computed 
by shifting the previous one, X 1 = X a + dX°, with 
dX° = Yli=i to obtain the new grid Gi - Now, 
the projector can be updated on Gi, and the MIG 
procedure applied for a second refinement. Itera- 
tions are terminated when, at any node, the norm 
of the defect of invariance | A| = \f — P (f) \ is suf- 
ficiently small in comparison to the norm of the vec- 
tor field |/|. 

Finally, it is important to discuss the projector 
P appearing in Eq. (fTU]) . The MIG method makes 
use of the thermodynamic projector [3], whose con- 
struction is briefly reviewed below. Let VG and r 
be the gradient of G and the tangent hyperplanc, 
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evaluated at a given grid node X 7 respectively. Let 
To = -rnfcer(VG), where fcer(VG) indicates the hy- 
perplane orthogonal to VG. Assuming that t ^ t , 
let U\ be a vector of the tangent plane r, such that 
1 and 



Miifcc 7 =0, H = 



8 2 G 



(11) 



_d[x i \d[x j ] 

where x is an arbitrary vector of the subspace To- 
The thermodynamic projector acts on a generic vec- 
tor r] as follows 



Pr] = (t7VG t ) ui 



En 
i=2 



Here, the set of vectors {1x2, 
of To, such that 

• « T 



V Huj j fi, (12) 
, u n } forms a basis 



u i Hu i = Sij, Vi,j = 2 



, n 



(13) 



with Sij denoting the Kronecker delta. In the case 
t = to, let {mi, . . . ,u n } be a basis of r such that 
UiHuJ = Sij, then (fT2"]) takes the form: 



(14) 



It is worth noting here a remarkable feature of the 
thermodynamic projector: The construction of (|12p 
on the SIM performs a slow-fast motion decomposi- 
tion. In other words, close to the SIM, the slow dy- 
namics of the system takes place in the image of 
P, while the fast dynamics evolves in its null space. 
More details about the Jacobian matrix J, the gra- 
dient VG and the second derivative matrix H are 
discussed in the appendix A. 

3.2. Initial approximation of the SIM 



the solution of (fTS"]) is unique if it exists [5]. For 
our purpose, it proves convenient to consider the 
whole g-dimensional manifold of constrained min- 
ima, regarding the fixed quantities as parameters. 
Here, it is worth to point out a connection between 
the notion of QEM and the method of Rate Con- 
trolled Constrained Equilibrium (RCCE) [9J. The 
RCCE method assumes that the reduced dynamics 
of the system ([3]) evolves along the QEM, obtained 
with a special choice of the vector set I 3 . Typically, 
one constraint concerns the total number of moles, 
Sr=i t-^"*] = £ X ' wm l e others may be related, for in- 
stance, to active valence and free oxygen |10j . 

For our purposes, the QEM must be explicitly 
constructed in the concentration space and refined 
via the MIG procedure. Details about the special 
choice of the set of constraints employed here are 
given in the subsequent sections. The grid-based 
approximation of the QEM is adopted as the ini- 
tial grid for the iterative procedure MIG. This 
approximation is constructed following the quasi- 
cquilibrium grid algorithm (QEGA) [5] and pro- 
ceeds as follows: Starting from an initial point X° 
close to the QEM (e.g. the equilibrium itself), the 
function G is approximated by a second-order poly- 
nomial around X°, and the minimization problem 
(fT5|) is recast to an algebraic system. The solution 
of this system delivers a new node X 1 in the neigh- 
borhood of the former one and close to the QEM. 
Similarly, starting from X 1 , the procedure can be 
applied again for seeking further grid nodes, and 
it is terminated when negative concentrations are 
obtained. 



As suggested in [3J, the notion of a quasi- 
equilibrium manifold (QEM) can be used for initial- 
izing the MIG procedure (see also [5]). In general, 
a q-dimensional QEM represents a manifold in the 
concentration space which is given by minimizing 
the Lyapunov function ^ under a set of q linear 
constraints expressing a re-parametrization of the 
original variables [Xi] in terms of some new vari- 
ables f . Let us consider the minimization problem: 

min G 

" (15) 
s.t. ^l{[X i }=¥, 1 ,/ 

i=l 

where I 3 = l 3 n ) is a set of q fixed vectors 

that will be specified below. Because of the convex- 
ity of G, for each fixed value of the quantities £■? , 



4. Application of MIG to reduction of a 
detailed mechanism 

In the sequel, a 7?2-air system reacting according 
to the nine-species, 21-step detailed mechanism of 
Li et al. [11], is studied. An adiabatic constant vol- 
ume reactor with 7?2-air mixture in stoichiometric 
proportions is considered, where the density and the 
mixture-averaged specific energy are chosen as p = 
4.58 kg/m 3 and e = 1.28 MJ/kg, respectively. The 
one- and two-dimensional slow invariant manifolds 
are described by constructing the pertinent invari- 
ant grids, which are utilized to integrate the reduced 
system. In the sequel, the concentration of species k 
is expressed in terms of specific mole number, (f>k — 
[Xk]/p. 



4 



4.1. Thermodynamic projector 

In the present section, we explicitly discuss the 
construction of the thermodynamic projector for 
1-D and 2-D grids, in a nine-dimensional concen- 
tration space. Let a generic 1-D grid Q be given as 
a collection of points in the concentration space. 
Assuming that a parameter £ is uniquely associ- 
ated to any grid point, the tangent vector u = 
(d{Xi]/d£, ...,d[A 9 ]/d£) can be approximated at 
any node via finite differences. Any nine-component 
vector r) can be projected onto Q as follows: 

P(ri) = — ^(VGr] T )u (16) 

For a 2-D grid, two parameters and £ 2 are asso- 
ciated with each point, so that two tangent vectors 
Mi = [d[Xi]/ dt; 1 ) and u 2 = (d[Xi}/d^) can be 
evaluated at any grid node. In order to construct the 
thermodynamic projector, it is convenient to intro- 
duce the 11 x 12 block matrix 







VG 
—I 



(17) 



with I denoting the identity matrix. Let us assume 
that A is a full rank matrix, and 1*2 is formed by the 
first nine components of a vector spanning the null 
space of A. According to the notations introduced 
in section [3TT1 let r indicate the tangent hyperplane 
spanned by U\ and u^. The intersection tq = r n 
/cer(VG) is one-dimensional and U2 is a basis of To, 
that is, the 2-D thermodynamic projector acts, on 
an arbitrary vector 77, as follows: 

1 , m\ 1 



p(n) = 



VGu\ 



(VGr] T ) «! 



" T " T 

U2HU2 



T]HU2 J U2 

(18) 

where the vector Ui is parallel to the hyperplane r, 
such that UiHu^ = 0. If A is not of full rank, we 
take U2 — u 2 and U\ parallel to r with UiHuT^ = 0. 



4.2. Construction of invariant grids 

For the case under study, a 1-D spectral quasi- 
cquilibrium grid (SQEG) was constructed as the 
first approximation of the 1-D SIM. In particular, 
the vector I 1 appearing in (|15p was taken as the 
left eigenvector of the Jacobian matrix J evalu- 
ated at the equilibrium point corresponding to the 
eigenvalue with the smallest absolute value [516] . 
The initial grid, as well as any subsequent grid, is 






























\ 




\ 



Fig. 1. Starting from the equilibrium point (filled circle), 
the 1-D SQEG (diamonds) was constructed via QEGA and 
refined via MIG to obtain the 1-D invariant grid (squares). 



Intermediate grid 




0.02 

0.01 <> [mot/kg] 



Fig. 2. The 2-D SQEG (continuous lines) is constructed and 
refined via MIG iterations. A projection of an intermediate 
refined grid (dots) is shown. 

parametrized by £ = £? =1 1\ [JQ]. The 1-D SQEG 
was refined via MIG till the dimensionless ratio 
between the norm of the defect of invariance and 
the vector field |A|/|/| became smaller than 0.001 
at every grid node. The results are shown in Fig. 
[TJ Here, it is worth to mention that the SQEG and 
the invariant grid are in a good agreement in the 
neighborhood of the equilibrium point. Moreover, 
the SQEG also proves to be a good approximation 
with respect to the major species in the full concen- 
tration space (as can be seen by its projection in 
the 4>h2-4>h 2 o subspace, Fig.Q})). 

A 2-D SQEG was also constructed by solving the 
minimization problem ()15j) . where the two vectors 
I 1 and I 2 were chosen as the two left eigenvectors 
of the Jacobian matrix evaluated at the equilibrium 
point corresponding to the two smallest eigenvalues 
in absolute value. The two reduced variables associ- 
ated with the grid nodes were ^ = J2^=i ^\ [Xi] and 

= ELi l i i X i\- The 2 " D S Q EG was a S ain refined 
until the threshold value 0.001 for the ratio | A\/\f\ 
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Fig. 3. The 2-D invariant grid (thin lines) was computed by 
refining the pertinent 2-D SQEG. The 1-D invariant grid 
(squares) and solution trajectories (bold lines) are reported. 

was reached. If the defect of invariance at a refined 
node kept increasing after several iterations, the new 
node was discarded. The 2-D SQEG accurately de- 
scribes the invariant grid only near the equilibrium 
point, and several Newton iterations (flO|) were re- 
quired during the refinement process. Figure[2]shows 
both the initial SQEG grid (solid lines) and its re- 
finement after three iterations. 

The projection of the final 2-D invariant grid onto 
the 4>h-4>H02-4>h 2 02 subspace is shown in Fig. [3l 
For the construction of both the 1-D and 2-D ther- 
modynamic projector, approximation of the tangent 
vector Uj = (<9 by first-order finite differ- 

ences was found to be sufficient. 

4.3. Reduced system 

Once the invariant grid is obtained, it can be 
stored in tables and used during the time integration 
of the reduced system. Indeed, let us assume that 
the 1-D invariant grid Qt nv is constructed and a vec- 
tor m — (mi, ...,mg) is chosen in such a way that 
the parameter £ = 5^»=i m i[^i] i s uniquely associ- 
ated with every point of Qinv . The original system 
Q reduces to the single equation: 

|=P(/(e))m r (19) 

When a 2-D reduced description is adopted, two vec- 
tors are introduced (m 1 , m 2 ) so that the new vari- 
ables are = Yn=i m \[ x i]^ 2 = Yn=i m il x i] and 
the reduced system reads: 

dt (20) 




0.5 1 1.5 2 2.5 3 




Fig. 4. Starting from a point located on the 2-D invariant 
grid, the reduced system was integrated by using an ex- 
plicit Runge-Kutta 4-th order scheme with a fixed time step 
At = 10 -8 s (symbols). Continuous lines represent the so- 
lution of the detailed model. 

Notice that the dependence of the vector field / 
on the reduced variables is not explicitly known and 
a proper look-up table is needed during time inte- 
gration of (JT9J) or ([20)) . In other words, a continua- 
tion procedure of the approximate invariant mani- 
fold from the discrete grid has to be implemented. 
For the problem under study, the reduced variables 
of the 2-D invariant grid were chosen according to 
the SQEG parameterization: m 1 = I , m 2 = I 2 , and 
the continuation of the invariant manifold from a 
generic 4-node cell of the invariant grid Qi nv was ob- 
tained by linear interpolation. Since the chosen pa- 
rameterization of Qi nv leads to a non-regular Carte- 
sian grid in the parameter space, the 4-node cell 
was mapped to a standard rectangle where a bi- 
variate linear interpolation was used to reconstruct 
the point on the SIM. The system (|20|) was solved by 
an explicit 4-th order Runge-Kutta scheme with the 
time step At — 10 -8 s. The results were compared 
with the solution of the detailed system ([3]), obtained 
with the same ODE solver. However, in the latter 
case the time step needed was one order of magni- 
tude smaller due to the stiffness of the detailed sys- 
tem. The comparison, shown in Fig. 2J proves both 
that (|2"01) is less stiff than © and a linear interpola- 
tion on the 2-D invariant grid in Fig. [3] delivers an 
excellent approximation of the correspondent slow 
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invariant manifold. The results were compared on 
the basis of the relative deviation of the reduced so- 
lution with respect to the detailed one averaged in 
time. The maximum error was found to be around 
2% for the evolution of 4>h 2 02 > while for the remain- 
ing species the mean relative deviation was below 
1%. 

Finally, it is worth to point out the computational 
effort needed for the construction of the invariant 
grid. For the case under study, the initial 2-D SQEG 
contains 1650 nodes. It was generated in 6 seconds 
and refined in about 10 minutes on a single processor 
3 GHz by using a Matlab code. 



5. Conclusions 

In this work, the Method of Invariant Grids (MIG) 
is applied for the first time to reduce a detailed 
hydrogen mechanism in non-isothermal conditions. 
The two-dimensional reduced model was then com- 
pared to the detailed one in an adiabatic constant 
volume reactor with H2-3at in stoichiometric pro- 
portions. 

The Spectral Quasi Equilibrium Grid (SQEG) [5] 
proves suitable for providing the MIG with an ini- 
tial collection of points (initial grid). The one- and 
two-dimensional SQEG grids describe quite well the 
dynamics of the major species, but they are not 
able to capture the correct evolution of some of 
the radicals. Therefore, several MIG iterations are 
needed in order to construct accurate 1-D and 2-D 
discrete approximation of the slow invariant mani- 
fold (SIM), providing the reduced description of the 
original nine-dimensional system. A bi-variate lin- 
ear interpolation was used to reconstruct the SIM 
from the invariant grid during the time integration 
of the reduced system. As exemplified by the re- 
duction in the number of time steps needed to in- 
tegrate the reduced system by an order of magni- 
tude, the stiffness is significantly reduced, and the 
two-dimensional models proves to be an excellent 
approximation of the detailed dynamics. 
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7. Appendix A 

The gradient V67, the second derivative matrix 
H and the Jacobian matrix J in section [3~T1 can be 
written more verbosely as follows: 



VG 

H = 
J = 



dG 

d 2 G 



d[Xi]d[X 3 ] 
dfi 



d[X s 



(21) 



where the notation indicates that the partial deriva- 
tives are computed under fixed e and species concen- 
trations. The Jacobian matrix J acts on a generic 
vector rj as follows: 

r 
s=l 

(22) 

At the equilibrium point, the matrix J reduces to 

J '&3) = ~E + n r)7.(<) (Htf) u). 

s=l 

(23) 

The latter operator is symmetric in the following 
sense: 

r]J'Hv T = vJ'Hrf, (24) 
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where r\ and v are two arbitrary n-component vec- 
tors. Because of the symmetry of J', by using the 
latter matrix in Eq. (|10[) at any node, instead of the 
full Jacobi matrix J, the stability of MIG iterations 
can be improved [7] . 

Assuming that d elements participate in the reac- 
tion, and Nk is the total number of atoms k in the 
system, let ([Zi], [Z n -d-i]) be an independent 
subset of ([Xi], [X„]), so that the latter variables 
are linearly dependent on the former ones: 

[Xi] = [X t ]{[Z 1 ],...,[Z n _ d _ 1 ],p,N 1 ,...,N d ). 

The MIG refinements can be carried out directly in 
the subspace described by [Zi] . The Eqs. (fT0|) remain 
valid with the Lyapunov function in the form ([5]). 
Now, (|21|) are substituted with 

' dG \ 



VG = 

H 
J 



d[Zi]j 
d 2 G 



e,p,Ni,...,N d ,[Z jlH \ 



d[Zi]d[Z 3 ] 

dfi 
d[Z L 



(25) 



e,p,Ni,...,N d ,[Z k ^] 

and the set of vectors {bi, b^}, appearing in dTUJ) , 
represents a basis in the null space of the thermody- 
namic projector (|12|) . 

8. Appendix B 

In this work, calculations were carried out by us- 
ing the reaction mechanism suggested in [11] . here 
reported in Table [1] 



Table 1 

Detailed i?2-air reaction mechanism. Units are cm 3 — mol — 
sec - Real - K, and k 3 (T) = A s T Us exp(—E 3 /RT). a Troe 
parameter is: F c = 0.8. Efficiency factors are: £h 2 = 10.0, 
eh 2 = 1-0 and eo 2 = —0.22. ^Troe parameter is: F c = 0.5. 
Efficiency factors are: £h 2 = H-0, Eh = 1-5- 



Reaction 


As 


n s 


Es 


1. H 2 +0 2 ^O + OH 


3.55 x 10 15 


-0.41 


16.6 


2. O + H 2 ^ H + OH 


5.08 x f0 4 


2.67 


6.29 


3. H2+OH ^ H 2 + H 


2.16 x 10 8 


1.51 


3.43 


4. + H 2 O^OH + OH 


2.97 x 10 6 


2.02 


13.4 


5. H 2 +M^±H + H + M 


4.58 x 10 19 


-1.40 


104.38 


6. O + O + M ^± 2 + M 


6.16 x 10 15 


-0.50 


0.00 


7. O + H + M ^±OH + M 


4.71 x 10 1S 


-1.0 


0.00 


8. H + OH + M ^ H 2 + M 


3.8 x 10 22 


-2.00 


0.00 


9. H + O a (+M) ^ H0 2 (+M) a 


k 6.37 x 10 20 


-1.72 


0.52 


10. HO2 + H ^±H 2 + 2 


1.66 x 10 13 


0.00 


0.82 


11. H0 2 + H ^±OH + OH 


7.08 x 10 13 


0.00 


0.30 


12. H0 2 + ^0 2 + OH 


3.25 x 10 13 


0.00 


0.00 


13. HO2 + OH ^ H 2 + O2 


2.89 x 10 13 


0.00 


-0.50 


14. HO2 + HO2 ^H 2 2 + 2 


4.20 x 10 14 


0.00 


11.98 


15. HO2 + HO2 ^ H2O2 + O2 


1.30 x 10 11 


0.00 


-1.63 


16. H 2 2 (+M) ^±20H(+M) b 


k Q 1.20 x 10 17 


0.00 


45.5 




hoc 2.95 x 10 14 


0.00 


48.4 


17. H2O2 + # ^ i? 2 + OH 


2.41 x 10 13 


0.00 


3.97 


18. H2O2 + H ^ HO2 + H 2 


4.82 x 10 13 


0.00 


7.95 


19. H2O2 + ^OH + HO2 


9.55 x 10 6 


2.00 


3.97 


20. H2O2 + OH ^ HO2 + H2O 


1.00 x 10 12 


0.00 


0.00 


21. H2O2 + OH ^ HO2 + H 2 


5.8 x 10 14 


0.00 


9.56 
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